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We develop a new method for taking into account the interference contributions to proton-proton inelastic 
cross-section within the framework of the simplest multi-peripheral model based on the self-interacting scalar 
(j) 3 field theory, using Laplaces method for calculation of each interference contribution. 

We do not know any works that adopted the interference contributions for inelastic processes. This is 
due to the generally adopted assumption that the main contribution to the integrals expressing the cross 
section makes multi-Regge domains with its characteristic strong ordering of secondary particles by rapidity. 
However, in this work, w e find what kind of space domains makes a major contribution to the integral and 
these space domains are not multi-Regge. We demonstrated that because these interference contributions are 
significant, so they cannot be limited by a small part of them. With the help of the approximate replacement 
the sum of a huge number of these contributions by the integral were calculated partial cross sections for such 
numbers of secondary particles for which direct calculation would be impossible. 

The offered model qualitative agrees with experimental dependence of total scattering cross-section on 
energy y/s with a characteristic minimum in the range ^/s ^10 GeV. However, quantitative agreement was 
not achieved; we assume that due to the fact that we have examined the simplest diagrams of <p 3 theory. 

Keywords: inelastic scattering cross-section, total scattering cross-section, Laplace method, virtuality, multi- 
peripheral model, Regge theory 



I. INTRODUCTION 

This paper is the sequel to [1 and 2], where to calculate 
proton-proton scattering partial cross-sections within the 
framework of multi-peripheral model the Laplace method 
was applied. 

The inelastic scattering amplitude with production of 
a specified multiplicity of secondary particles, in frame- 
work of the multi-peripheral model can be represented 
as a sum of diagrams demonstrated on Fig.l. To calcu- 
late the partial cross-section <j n is necessary to evaluate 
an integral of the squared modulus of a sum of contribu- 
tions shown in Fig.l. After simple transformations [2], 
the expression for the partial cross-section can be rep- 
resented as a sum of "cut" diagrams in Fig. 2. We call 
summands entering into the sum Fig. 2 the interference 
contributions. Approximate calculation of their sum is 
the purpose of this paper. 

At present time the inelastic scattering processes are 
considered without the interference contributions [3 and 
4]. This due to the generally adopted assumption that 
the main contribution to the integrals expressing an in- 
elastic processes makes multi-Regge domains [3-6] with 
its characteristic strong ordering of secondary particles 
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FIG. 1. Diagram representation of an inelastic scattering 
amplitude when the n secondary particles are formed. Here 
Pi and P2 are the four-momenta of primary particles before 
scattering; P3 and P4 are the four-momenta of primary parti- 
cles after scattering; Pi 1 ,Pi 2 , • • • ,Pi n are the four-momenta of 
secondary particles. Symbol ^ denote a sum over 

P(*l, 12,---, in) 

all permutations of indices i\ — 1, 12 =2,...,i n — n. 



by rapidity. This means that the rapidity of neighboring 
particles on the "comb" should be different from each 
other by a large value. Thus the amplitude of the right- 
hand and left-hand parts of the diagram on Fig. 2 for 
different orders of connecting lines would be significantly 
different from zero to almost non-overlapping regions of 
phase space and integral of their product would be a 
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FIG. 2. Representation of the partial cross-section as a sum 
of "cut" diagrams. The order of joining of lines with four- 
momenta pk from the left-hand side of the cut is as following: 
the line with p\ is joined to the first vertex, the lines with p2 is 
joined to the second vertex, etc. The order of joining of lines 
from the right side of cut corresponds to one of the n\ pos- 
sible permutations of the set of numbers 1, 2, . . . , n. Where 
Pj(k), k = 1, 2, . . . , n denote the number into which a number 
k goes due to permutation Pj. An integration is performed 
over the four-momenta pk for all "cut lines" taking into ac- 
count the energy-momentum conservation law and mass shell 
condition for each of pk ■ 



small quantity. 

However, as it was shown in [1] near the threshold of 
the n particles production at the maximum point of the 
scattering amplitude Fig.l difference between neighbor- 
ing particle's of rapidities is close to zero and at higher 
energies increases logarithmically with energy y/s growth. 
This difference has factor l/(n+l), so for high numbers of 
secondary particles it increases slowly with energy. More- 
over, even if each of interference terms is insignificant, all 
of them are positive and a huge amount n! of them not 
only makes it impossible to discard them, but also leads 
to the conclusion that the contribution of a "ladder" di- 
agram Fig. 2, which is usually only taken into account, is 
negligibly small compared with the sum of the remaining 
interference terms. This was shown in [2]. For the rel- 
atively small number of secondary particles (n < 8) we 
are able to calculate all the interference contributions in 
the direct way without any approximations. 

Further in this paper we will demonstrate method for 
approximate calculation of the sum of the interference 
contributions for large numbers of secondary particles, 
when direct numerical calculation is not feasible. 



II. METHOD DESCRIPTION 

Using the Laplace's method we have found [1 and 2] 
the mechanism of partial cross-section growth, which was 
not taken into account in the previously known vari- 
ants of multi-peripheral model. This mechanism may 
be responsible for the experimentally observed increase 
of hadron-hadron total cross-section. However, in this 
approach based on the Laplace's method, it was found 



out that the calculation of partial cross-sections in the 
multi-peripheral model can be limited just to contribu- 
tions from the "cut ladder diagram". Because for any 
number of the secondary particles n there is the wide 
range of energies y/s, where such contribution is negligi- 
bly small compared to the sum of n\ positive interference 
contributions. At the same time, as we will demonstrated 
further, the allowance for the interference contributions 
results in the appearance of multipliers in expression for 
the partial cross-section, which are decrease with the en- 
ergy rise (see below Eq.7). Thereupon the question 
arises: "Will the sum of partial cross-sections increase 
with energy rise if we take interference summands into 
account?" 

As shown in [1], each term in sum shown in Fig.l with 
accuracy up to the fixed factor is a function with real and 
positive values, which has a constrained maximum if its 
arguments satisfy the mass-shell conditions and energy- 
momentum conservation law. Therefore, in the c.m.s. of 
initial particles function corresponding to the left-hand 
part of cut diagram in Fig. 2 can be rewritten in the neigh- 
borhood of maximum point in the form [1 and 2] 



a(x)=a (jf (°>) 



where X is the column composed of 3n + 2 indepen- 
dent variables, on which the scattering amplitude de- 
pends after consideration of mass-shell conditions and 
energy-momentum conservation law; the first n compo- 
nents of column are the rapidities of secondary particles; 
the next n components are the x components of transver- 
sal momenta of secondary particles (it is supposed that 
the reference system is chosen so that Z-axis is directed 
in the line of the three-dimensional momentum Pi of ini- 
tial particle in Fig.l), the y - components of secondary 
particle transversal momenta and the two last variables 
are the antisymmetric combinations of particle transver- 
sal momenta P 3 and P 4 , i.e., X 3n+1 = \ (P 3 ± x ~ Pa±x) 
-^3n+2 = \ (P3±y — P^±y) We denote the column of the 
values of variables in a maximum point through and 
a matrix with the elements 



D ab 



d 2 

dX a dX t 



(„(,(*))) 



X=X(°) 



(2) 



where 



a = 1, 2, • • • , 3n + 2, b = 1, 2, . . . , 3n + 2 (3) 

are the coefficients of the Taylor series expansion of 
amplitude logarithm in the neighborhood of maximum 
point. As it was shown in [1], if we do our computations 
in the c.m.s. of initial particles, the maximum is reached 
when transversal momenta is zero and secondary particle 
rapidities are close to numbers that formed an arithmetic 
progression. 

If we denote the difference of this progression through 
Ay (n, y/s) and the value of particle's rapidity to which 
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the line attached to the fc-th vertex of diagram in Fig.l 
corresponds, through 

Ay(n,^) = yW - y^, k = 1, 2, • • • , n - 1 (4) 
we get [1]: 



Vk 



(0) 



(s±l - k ) Ay(n,y/S), k = l,2,- 



(5) 



The form of the function Ay (n, y^) has been discussed 
in [1]. For further consideration, it is important that it is 
a slowly increasing function on s and decreasing function 
on the number n of the secondary particles and vanishes 
when s is equal to the threshold of n particle production. 
Thus, the column contains only the first nonzero n 
rapidity components, which are defined by Eq.5. 

The following expression corresponds to the right-hand 
part of cut diagram in Fig. 2: 

P j (a(x))=a(x^) 

xexp (PjX - X^D (PjX - X^^j (6) 

The interference contribution corresponding to whole 
"cut" diagram, which correlates with the j-th summand 
in Fig. 2, is proportional to an integral of the product of 
functions Eq.l and Eq.6 over all variables. Denoting an 
interference summand corresponding to the permutation 

Pj through (j' n {P^j an d calculating its Gaussian integral 

(at the same time, other multipliers besides the squared 
modulus of scattering amplitude in an integrand are ap- 
proximately replaced by their values at the maximum 
point [2]), we get 



yj det (i (p 
xexp(-i((Axf ) T J D«AX 
where we use the following notations: 



PJDP 3 



(0) 



(7) 



AX) 



(0) 



X u -P~ 



i (1(0)) 
v(y/a) = l — 



(8b) 
(8c) 



E 



^~s^s/i-M^{E P /2) ^(E P /2) 2 -M 
n 

P = vs-J2 ch {yi 0) ) ( 8d ) 



k=i 



M is the mass of initial particle, which is made dimen- 
sionless by the mass of secondary particle (it is supposed 
that the energy y/s is also made dimensionless by the 
mass of the secondary particle). 

Note, that here and in the following sections we will 
use the "prime" sign in ours notation to indicate that 
we use a dimensionless quantity that characterized the 



dependence of the cross-sections on energy, but not their 
absolute values. 

The value of amplitude at the maximum point 
increases with the yfs growth due to mecha- 
nism of virtuality reduction [1]. However, the distance 
AX^ between maximum points of "cut" diagram also 
increases with the y/s growth. Therefore, the exponential 
factor entering in Eq.7 can decrease with energy growth. 
This makes considered above question. How competition 
of these two multipliers will result on the dependence of 
the sum of partial cross-sections on yfsl 

Thus, each interference contribution can be computed 
numerically. However due to the huge number of con- 
tributions and large number of secondary particles n the 
direct numerical calculation of the sum of interference 
terms in Fig. 2 is impossible. We can avoid this difficulty 
in the following way. The maximum in the right part 

of cut diagram in Fig. 2 is attained at X = P- 

-i (1(0)). 

In other words, a maximum of function, which is asso- 
ciated with the right-hand part of cut diagram, can be 
obtained from a maximum of function, which maps with 
the left-hand part of cut diagram, by the rearrangement 
of arguments. Then the value of each interference con- 
tribution is determined by the distance between points 
of maximum in the right-hand and left-hand part of cut 
diagram as well as by the relative position of these max- 
imum points, since in different directions contributions 
to scattering amplitude fall off with distance from point 
of maximum, in general, with different rate, and also by 
the relative position of proper directions of the matrices 
D and PjDPj. In other words, multiplying Gaussian 
functions corresponding to the right-hand and to the left- 
hand part of interference diagrams in Fig. 2 each time we 
will obtain as a result Gaussian function, which has the 
proper value at the maximum point (which we call the 
"height" of the maximum) and the proper multidimen- 
sional volume cutout by resulting Gaussian function from 
an integration domain (which we call the "width" of the 
maximum) . 

We assume that summands in Fig. 2 are arranged in 
ascending order of the distance between the maximum 
points in the right-hand part and left-hand part of cut 
diagram (we denote this distance through r) so that 
" cut" diagram with the initial attachment of lines to the 
right-hand part of diagram corresponds to j = 1. In 
other words, the line of secondary particle with the four- 
momentum pi is attached to the i-th top in the right-hand 
part of cut diagram in Fig. 2. As follows from Eq.7, the in- 
terference contributions exponentially decrease with the 
r 2 growth. However, in spite of this the interference con- 
tributions do not become negligible due to their huge 
number, which, as discussed below, are increases very 
rapidly with r 2 growth. The value of r 2 is proportional 
to the square of magnitude Ay(n,y/s), which, as was 
noted above, is zero on the threshold of n particle produc- 
tion and slowly increases with distance from this thresh- 
old. Therefore, for each number n there is the fairly 
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FIG. 3. The interference contributions dependence on zi at 
y/s = 1000 GeV: (a) n = 8, (b) n = 9. Here and in sub- 
sequent figures the interference contributions divided by the 

/ 3n+2 3n+2 

common multiplier exp I — ^ 

\ a=l b=l 

cated on the y-axis. Obviously, that to the one value of zi 
correspond a lot of different contributions, as well as that the 
average values of the logarithms of these contributions are 
placed approximately on a straight line (see below Eq.25 and 
Fig.4). 



E X ( a 0) D ab X ( b 0) ) are indi- 



wide range of energies close to the threshold, in which 
the sharpness of decrease of the interference contribu- 
tions with the r 2 increase is small in the sense that it is 
less important factor than the increase in their number. 
At such energies, which we call "low", the partial cross- 
section o' n is determined by the sum of huge number of 
small interference contributions. When the magnitude 
Ay(n, y^) is increased with the further growth of en- 
ergy y/s, the decrease rate of interference contributions 
increases, while the growth rate of their number with the 
r 2 increase does not change with energy. At such ener- 
gies, which we call "high" , the main contribution to the 
partial cross-section is made by the relatively small num- 
ber of interference terms corresponding to the small r 2 , 
which can be calculated by Eq.7. 

If we compose the n-dimensional vector (we denote it 
through y^ ^) from the particle rapidities Eq.5, which 



constrainedly maximizes the function associated with 
the diagram with the initial arrangement of momenta 
in Fig. 2, vectors maximizing the functions with another 
momentum arrangement will differ from the initial vec- 
tor only by the permutation of components, i.e., these 
vectors have the same length. Consider two such n- 
dimensional vectors, one of which corresponds to the ini- 
tial arrangement, and another - to some permutation, 
then in the n-dimensional space it is possible to "pull 
on" a two-dimensional plane on them (as a set of their 
various linear combinations), where two-dimensional ge- 
ometry takes place. Therefore, the distance r will be 
determined by cosine of an angle between the consid- 
ered equal on length n-dimensional rapidity vectors in 
the two-dimensional plane, "pulled" on them. An an- 
gle corresponding to the Pj permutation we designate 
through 0j, < 0j < tt. 

Thus, each of the terms in the sum Fig. 2 can be 
uniquely matched to its angle 0j. At the same time the 
variable z = cos (6) is more handy for consideration than 
an angle 0j. Using Eq.5, can be shown that the variable 
z can take discrete set of values: 



1 



12 



* = 0,1, 



(n- l)n(n + l) 
(n — 1) n (n - 



6 



(9) 



Note that although the relation Eq.5 for the rapidities 
of secondary particles is satisfied with high accuracy at 
the maximum point, it is still approximate. This means 
that those contributions, to which matched one and the 
same value of variable z in Eq.5, in fact, matched a 
slightly different from each other values of z. 

Consequently, to such contributions correspond a sim- 
ilar but unequal to each other distances between maxi- 
mum points in a "cut" diagram. In addition, this dis- 
tance, as was discussed above, is not a unique factor af- 
fecting to the value of interference contribution. There- 
fore, if each interference contribution is associated with 
the value of variable z by the approximation Eq.5, it 
appears, that the different values of interference contri- 
butions correspond to the one and the same value of z\ 
(see Fig. 3). 

Thus, while each contribution is associated to some 
value of variable z in the approximation Eq.5, the value 
of contribution is not the unique function of z. However, 
the sum expressing the partial cross-section cr' n can be 
written in the following way 



(n-l)n(n + l) 



( 



i=o 



AN, 



V 



(10) 



where ANi the number of summands to which the value 
Zj = z\ is corresponds in the approximation Eq.5. The 
average value of all interference contributions in Eq.10 is 
already the unique function of z\. Therefore, we intro- 
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FIG. 4. Comparing the values of In ((cr' n (zi))) obtained by a direct numerical calculation with consideration of all interference 
contributions (circles) and by approximation Eq.24 (straight line) at n = 8, = 10 GeV (a); n — 9, y/s — 10 GeV (b); n — 8, 
yfi = 100 GeV (c); n = 9, = 100 GeV (d). 



duce notation 



E < (Pi) 



= «(«)) 



(11) 



where {a' n (zi)) is some function, whose form at "low" 
energies can be determined from the following consider- 
ations. 

For any multiplicity n when the values of parameter 
I in Eq.9 are small and when number of corresponding 
interference contributions is relatively small, we can di- 
rectly calculate these elements and their sum. Denote 
the maximum value Z, for which all interference contri- 
butions are calculated through i - In particular, in this 
paper we managed to calculate the interference contribu- 
tions up to Iq = 6. Partial cross-section can be written 
as 

a > = a f(h) + /(0 = 

u n u n * u n 

(n-l)n(n + l) 

= E < (Pj) + t AiV ; (a' n ( Zl )) (12) 

«=0,1,— / 
Ah) 



energies, and affl is the sum of contributions sufficient at 
"low" energies. Thus, the difficulties in the calculations 
of the huge number of interference contributions mainly 
relates to the range of "low" energies and can be reduced 
to the approximate calculation of (<J f n (zi)) and AiVj. 



III. THE APPROXIMATE CALCULATION OF {a' n (zi)). 

As follows from Eq.7, the exponential factor exerts the 
most significant effect on the dependence of (a f n (zi)) on 

z t . Note that the expression (AXj 0) ) D^AX^ en- 
tering into the exponent in Eq.7 depends only on those 
matrix components, which are at the intersection 
of the first n rows and first n columns, since all column 
AXj ^ components starting with n + 1 are zero, because 
they are the particle momentum transverse components 
at the maximum point. If we denote the matrix com- 
posed of elements located at the intersection of the first 
n rows and first n columns of the matrix through 



where <r n is the sum of contributions sufficient at "high" Dy and a matrix, which is obtained from the matrix D 
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FIG. 5. The values of (w (zi)) obtained by direct calculation values Eq.24 for all interference contributions for n = 8 and n = 9 
at y/s = 10 GeV 5(a), 5(c) accordingly; for the same number n, but at y/s — 100 GeV 5(b), 5(d) and the ratio (w(zi) / (w(zi ) 
for n = 8 and n = 9 at yfe = 10 GeV 5(e),5(g); y/s = 100 GeV 5(f),5(h). 



in analogy, through D yi we have 

m : D« = (f)- 1 + PjD- l P 3 ) _1 (13) 

The matrices D~ x and P^D~ 1 Pj have one and the 
same eigenvalues, but they correspond to different eigen- 
vectors. We denote the normalized to unit eigenvec- 
tor corresponding to the minimal eigenvalue of matrix 
Dy 1 + PJ D~ 1 Pj through u m i n and the eigenvalue itself 
- through A m i n . This implies 

Amin = ^min y ^min ~l~ ^min j y ^j'^min (1^) 

Since the minimum eigenvalue of matrix Dy 1 is equal 
to the minimum values of quadratic form D~ x u for 
the unit vectors u, the magnitude u^^Dy 1 Uuim is not 
less than the minimum eigenvalue of matrix Dy 1 . By 
analogy the magnitude u^ in Pj D~ x PjU m i n is not less 
than the minimum eigenvalue of matrix PJ D~ 1 Pj, which 

coincides with the minimal eigenvalue of matrix Dy 1 
and is reciprocal of the maximum eigenvalue of matrix 
D y denoted through d™ ax . Thus, A min > From 

this it follows that, the maximum eigenvalue of matrix 
{t)- 1 + PfD^P^j does not exceed d™ ax /2. By anal- 
ogy we obtain that the minimum eigenvalue of matrix 



[by 1 + PjD^P^j 1 is no smaller than d™ in /2, where 
<i™ is the minimum eigenvalue of matrix D y . Thus, an 

interval enclosing the eigenvalues of matrix D y J is, at 
least, twice smaller than an interval enclosing the eigen- 
values of matrix D y . We can demonstrate that at ap- 
proximation of an equal denominators [1] the value of 
d™ ax can be estimated in the following way 




i.e., an interval enclosing the eigenvalues of matrix D y ^ 
at any energies and number of particles is less than unity, 
whereas at the considerable values of Ay(n,s), i.e. at a 
distance from the threshold, this interval is much less 
than unity. 

Therefore, if we reduce matrix D y ^ to diagonal form, 
it will be close to a matrix multiple of unit matrix. If we 
represent this matrix in the form 

= ^Sp (i)^ E + AbW (16) 

where E is unit matrix, the eigenvalues of the traceless 
matrix Abffl will be small. Then 
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2V Axf) D^Axf 



(l-cos(%)) 



2 ^ 

fe=i 



Ad 



,(°) 



(17) 



where Ad^'jj. are the eigenvalues of matrix AD y ^ , Vk n is 
the transformation matrix to the basis composed from 
the eigenvectors of matrix AD y ^ (the summation over 
reheated indices is supposed). The second term in this 
sum is small in comparison with the first one due to the 
smallness of eigenvalues Ad y ^\ as well as due to their dif- 
ferent signs (since the trace of matrix AD y ^ is zero, the 
different terms over k partially compensate each other). 
Therefore, we can adopt the following approximation: 



AX 



(o) 



ksp (A 



|y< 0) | 2 (i 



cos(%)) 



(18) 



To approximately calculate the trace of matrix D y ^ 



we select the spherically symmetric part of matrix D y 
representing it in the form 



±Sp(D v ) 



E + AD V 



(19) 



The results of numeral calculation of the eigenvalues 
of matrix D y (which are denoted through , k = 
1, 2, • • • , n) are shown in Table. I. It is obvious that most 
eigenvalues are close between themselves with the ex- 
ception of a few eigenvalues, which are substantially 
smaller. Therefore, these smallest eigenvalues have the 
highest absolute value of deviations from mean eigen- 
value ^Sp (^Dy^j. Since all the eigenvalues of matrix D y 

are positive, the deviation of eigenvalues from average 
value is less than this average in absolute value (see Ta- 



bid). Note that the matrix D y 
the following form: 



can be represented in 



£>P = ±s p (d v ) [e 




ADy + PjADyP, 



\Sp(by) 



E- 



(20) 



By analogy we can conclude that the minimum eigenvalue 
of matrix 



Aby+PjAbyPj 



(21) 



(which is maximum in absolute value, see Table. I) is 
greater than the doubled minimum eigenvalue of ma- 
trix AD y . This means that all the eigenvalues of matrix 

ADy+PfADyPj 



\S P (Dy) 



are less than unity in absolute value. It 



applies equally to the eigenvalues of matrices — 



Sp(Dy) 



and P 1 i \ Pj • Therefore, we can represent the ma- 
trix Dy^ as the expansion in powers of TVTrTT- Since 
is traceless by definition, then a nonzero 



matrix 



contribution to Sp (jjy^ in addition to the term of 



"zero" order 



^S P (by) 



E can give terms starting with 



the second-order. As it follows from Table. I, the maxi- 
mum in absolute value eigenvalue of matrix — 



AD V 



creases with the energy growth. Therefore, we can expect 
that at "low" energies higher-order terms will make neg- 
ligibly small contributions. In such an approximation we 



have: 



(22) 



Let Eq.7 is taken in place of Eq.ll in approximation 
Eq.22, then we have 



{a' n { Zl ))=(A (x(°))) 2 u( v / s) 

\yW\ 2 Sp(D y ) 



x exp — 



2n 



(1 - zC, 



ANl z^ Zl ^(D+P-DP,)) 
Let us introduce the following notation 



(w 



(«)> = AkE 



Zj=Z t 



(23) 



(24) 



If we assume that multiplier (w (zi)} is weakly dependent 
on z\, we obtain 



y 



2 Sp(l)y) 



2n 



(25) 
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where z\ is the minimum value of z\ for which can 
be numerically calculated all interference contributions. 
Therefore, the magnitude (<j' n (zi )} can be directly cal- 
culated numerically. The results of numerical calculation 
of In ((a f n (zi))) over all interference contributions in com- 
parison with the results obtained by Eq.24 are demon- 
strated on Fig. 4, it follows that such an approximation 
is acceptable at "low" energies. 

Results shown in Fig. 4 confirm also our assumption 
that (w (zi)) weakly depends on z\. To analyze this de- 
pendence we turn to Fig. 5. It is obvious, that the mag- 
nitude (w (zi)} takes small values at "low" energies. This 
means that 



det 



\(b + pjbP 3 )) 



(26) 



takes large values at the same energies. Indeed, as it 
follows from the expression for the matrix l), Eq.26 tends 
to infinity on the threshold of n particle production, and 
this means that at threshold the phase space of physical 
area of the inelastic process with n particles production 
takes place is equal to zero. 

Because of symmetry with respect to direction inver- 
sion in a plane of transversal momenta the mixed second 
derivatives with respect to rapidities and transversal mo- 
mentum components are zeros. As a consequence, the 
determinant Eq.26 is equal to the product of the three 
determinants, first of which is composed from second 
derivatives with respect to rapidities, the second is com- 
posed from the second derivatives with respect to the 
transversal momentum ^-components and the third one 
is composed from derivatives with respect to the transver- 
sal momentum components. All the three factors tend 
to infinity at the threshold energy. As it follows from 
a numerical calculation, a matrix determinant composed 
from the second derivatives with respect to rapidities re- 
duced quite rapidly with energy growth. Matrix determi- 
nants composed from the second derivatives with respect 
to transversal momentum components also reduced, but 
in a wide energy range, they remain quite large. There- 
fore, the value of Eq.26 is great at all j. Since the function 
\j\fx varies slightly at the great values of argument, the 
function (w (zi)) weakly depends on z\. 

To estimate roughly the function (w (zi)) we can re- 
place it by the Taylor expansion taking into account just 
linear contributions. The expansion coefficients are found 
by the calculating of (w (zi)} for z\ close to 1 and (—1). 
In these cases the values of 



1/Wdet (l/2 (t) 



PJDP, 



(27) 



were obtained directly for all proper interference contri- 
butions, and after that we obtain the values of (w (zi)} 
by averaging using Eq.24. 



The values in Fig. 5 have been obtained by the direct 
calculation of 



A k 5 l/Jdet(l/2 (p + ffDPj 



(28) 




FIG. 6. A sphere S2 and figure F^\ (is shown by points). 
Basis in the four-dimensional space is chosen so that the one 
of vectors coincides with the vector e4=(|,|,|,|), and the 
three basis vectors of three-dimensional subspace, into which 
depicted sphere is embedded, are perpendicular to e^. 




FIG. 7. The partition of sphere S2 by shortest arcs joining the 
points of figure F±\ into the two "hexagonal" and one "tetrag- 
onal" regions 7(a); (b) areas, which is located on the borders 
of 4 or 6 points belonging to figure F±\ can be divided be- 
tween those points into figures of equal area; (c) whole sphere 
S2 is divided into figures of equal area, each of which contains 
the one point of figure F±\ one of these shapes are painted in 
white. 



with consideration of all interference contributions at dif- 
ferent y/s. 

So, we have the following expression instead of Eq.25 
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( CT n ( z l)) = ( a 'n ( z lo)) ( w + W 1 (1- Zi)) X exp 



U0)| 



2 5p (4) 



2n 



0/ -ZJ ) 



(29) 
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Therefore, considering vectors P- as the el- 

ements of n-dimensional euclidean space, which we de- 
note through E n , then the ends of all vectors P- 

are lie on the unit sphere embedded into the (n — 1)- 
dimensional subspace of E n . We denote this sphere 
through S n -2 and shape formed by the set of points in 

which the ends of vectors Pr 1 (j = 1,2, ... ,n!) 

come, denote through F n \. In particular, when n = 4 the 
sphere S2 and figure graphically look like in Fig. 6. 



FIG. 8. Diagrams, which correspond (n — 1) vectors 
P- 1 (y (0) ) closest to vector y (0) . 



where the coefficients wo and w\ are found by above men- 
tioned method. 



IV. APPROXIMATE CALCULATION OF THE AN t 
VALUES 

Let us turn to the new variables 



(o) 



a / /~\ / (n+l)n(n— 1) 

Ay (n, Vs) - ' — L 



(30) 
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where are determined by Eq.5, Y^ \ k = 1, 2, • • • , n 
are considered as the components of vector Y^°\ which, 
as it follows from Eq.30 is of unit length. 

Thus, the angle Oj between the vector yt°^ = 

(jJi*\ • • • , Vn^ and vector PJ 1 (y^ ^) obtained 
by the permutation of corresponding components is 
the same as the angle between the vector Y^ = 
(y{ 0) , F 2 (0) , • • • , F n (0) ) and vector pr 1 (y(°)) . More- 
over, as it follows from Eq.5 



J°) _ AO) (0) _ (0) m m m (0) _ _ (0) . 

k = 1, 2, • • • , n 



(31) 



It follows that all vectors P- are orthogonal to 

vector 



l/y/n,l/y/n,...,l/y/n 



(32) 



V 



n components 



We examine some geometrical properties of figure F n i 
at arbitrary n . If we apply the permutation transfor- 
mation component to all vectors in the n-dimensional 
space, where the vectors Y^ are primordially defined, 
the examined (n — l)-dimensional subspace as well as a 
sphere 5Vi-2 and figure F n \ go into themselves. As it fol- 
lows from the group properties of permutation group, the 
each point of figure F n \ can be obtained from any other 
point by some transformation Pj~ • This means that the 
configuration of the points of figure F n \ relative to each 
of these points must be identical, that can be clearly seen 
in Fig. 7(a). 



As it follows from Eq.31, besides the end of each vector 
p-i a figure F n \ contains also the end of vector 

^— PJ 1 (y(°^^J, i.e., a figure F n \ has a center of sym- 
metry, which coincides with the center of sphere S n -2- 
In this case, if we using point of F n \ form path from the 

point P- 1 (y(°)) to the point (^ (0) ))> then [t 

will be simultaneously formed a centro-symmetrical path, 
that leads from (-P^ 1 (^ (0) )) to Pf 1 (^ (0) ) of figure 
F n \. 



Joining these paths we will obtain the closed path, 
which "girdles" the sphere SVi-2- If we assume that there 
is such a "girdling" path, inside of which are concentrated 
all points of figure F n \, we would find that the figure F n \ 
has a "boundary" and "internal" points, that would con- 
tradict the fact that spacing of all points relative to each 
point of the F n \ should be the same. In other words, the 
points of figure F n \ must "crawl away" all over the sphere 
S n -2 and can not be concentrated on some area of the 
sphere. 



If we consider a vector Y^\ then closest to it are 
the vectors corresponding to permutations P^ 1 , Z = 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

z z z 



(a) (b) (c) 



FIG. 9. Comparison of the interference contribution distribution by the variable z — cos (B) (histogram) and the plot of 
function p(z) = dN ^ z g+ dz ) (solid line) at n = 8, Az = 0.1 (a); n = 9, Az = 0.1 (b); n = 9, Az = 0.05 (c). Here AN is the 
number of interference contributions corresponding to value of z in the proper interval of Az width. 




FIG. 10. Comparison of the values of right-hand side and left-hand side of approximate equality Eq.37 at n — 8 (a, b) and 
n = 9 (c, d). Circles are the values of ANi calculated with consideration of for all interference contributions; crosses are the 
values of function p (zi-i) Az from Eq.35. 



1, 2, ■ • • — l defined by the following relation The type of "cut" diagrams corresponding to such per- 
mutations is shown on Fig. 8. 

y(°) if k < I 

fe (0)'. ' At the same time, all the components of vector 

5 f 1 (^ (0) )) = "I Y %) ' ^ k = ^ 33 ^ (^ (0) ) ~ ^ (0) ' 6XC6pt /_th 1 + 1: ^ Z6r °' 

^ fc — Z + 1, whereas these two components take on the least values 

y fc (0) ,if k>i + i. 
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in modulus 



and 



( \l (n+l)n(n-l))' 



respec- 



(n+l)n(n-l) 

tively. 

Thus, we can conclude that the each point of figure 
F n \ has (n — 1) nearest neighboring points, which lying 
at distance of from it: 



24 



+ l)n(n- 1) 



(34) 



Connecting the each point of figure F n \ with its (n-1) 
nearest neighbors points by shortest arc thereby we di- 
vide the sphere 5^-2 into closed regions as is shown in 
Fig. 7(a). Indeed, let us choose the some point Aq of fig- 
ure F n \, and will move from it to the nearest point A\ 
along a shortest arc, then we move from the point A\ to 
the nearest point A2 etc. At the same time, motion in a 
backward direction is prohibited. Thus, there are (n — 1) 
paths going out from each point, and (n — 2) paths are 
allowed at each step. But since figure F n \ has the finite 
number of points at some step we will surely come back 
to the point Aq. 

Moreover, since shortest arcs joining two nearest points 
are subtended by equal chords r n in length (see Eq.34), 
this arcs are of the same length. Let us consider any 
two neighboring points Ai and A i+ i of figure F n \. Un- 
der any transformation PJ 1 the shortest arc, which joins 
the points Ai and and an arc joining the points 

PJ 1 (Ai) and PJ 1 (A i+ i) are of the same length. This 
means that the boundaries of closed regions formed by 
shortest arcs, which join neighboring points, replaced 
into one another under any transformation Pj 1 - It fol- 
lows that, if we examine closed areas which include any 
point of figure F n \, then the adjacent areas to all points of 
this figure will have the same "area" . There is one more 
requirement, to which the areas obtained by partition of 
the sphere 5Vi-2 must satisfy: they must not overlap, i.e., 
these regions do not have common internal points. In- 
deed, otherwise, at least any two of the examined arcs 
would intersect in some internal point of these arcs. As 
it follows from Eq.34, when n is large the value of r n 
is small. This means that when we join the each point 
of figure F n \ with its nearest neighbors by the shortest 
arcs of sphere 5Vi-2, these arcs practically coincide with 
chords, which tights them. 

If we assume, that any two chords A^A^+i and 
Ai 2 Ai 2 +i intersect in an internal point, then it is pos- 
sible "to pull" on them a two-dimensional plane. Then 
we get a flat rectangle A^A^ A^ 1+ i^ 2+ i, which has at 
least one angle no smaller than 90°. This means that 
square of diagonal lying opposite it is not less than sum 
of squares of the parties that make up the corner. Denot- 
ing the lengths of these sides through a and 6, we have 
a 2 + b 2 < r 2 . In this case, either a or b would not ex- 
ceed r n /\/2, i.e., the figure F n \ contains points, which are 
at distance less then r n but that cannot happen due to 
minimality of this distance. 

Thus, we can conclude that at an arbitrary n a sphere 
S n -2 can be divided into the parts of equal area, each of 



which contains only one point of figure F n \, as it shown 
in Fig.7(b)-7(c). 

Let us introduce a multidimensional spherical coordi- 
nate system so that the end of vector is the "north 
pole" of sphere S n -2. Then the number of points of fig- 
ure F n \, to which the values of variable z = cos (G) in 
the interval [z, z + dz] correspond, is equal 



dN (z, dz) = p (z) dz 



where 



(35) 



(36) 



T is the Euler gamma function. 

To verify the validity of Eq.35 we can calculate all in- 
terference contributions and corresponding values of z at 
n = 8 and n = 9 (since for the larger number of particles 
this can not be realized). 

The distributions of interference contribution from the 
variable z = cos (0) and the graphs of function p (z) = 
dN{z ' d z z +dz) from Eq.35 are shown in Fig.9. Obtained re- 
sults of numerical calculation of interference contribu- 
tions and by Eq.35 are in a good agreement. 

Moreover, as it follows from Fig. 9(b) and from Fig. 9(c) 
this fitness is improved with increasing number of parti- 
cles n, i.e., Eq.35 is suitable for large n, when the direct 
numerical calculation of all interference contributions is 
impossible. 

Taking Eq,35 and Eq.9 into account we obtain the fol- 
lowing the approximate equality 



where 



Az 
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(n-l)n(n-hl) 



(37) 



(38) 



Verification results of Eq.37 at n = 8 and n = 9 are 
presented in Fig. 10. 

Another verification of considered above equations is 

presented in Fig. 11, where the values of ^ a' n ( Pj J and 

Zj=Zl ^ ' 

approximating magnitudes (a' n (zi)) p (zi-\) Az (here 
(a' n (zi)) is calculated by Eq.29) are compared. 

From results demonstrated on Fig. 4 and Fig. 10-12, we 
can conclude that the at least for those numbers of parti- 
cles for which it can be directly tested Eq.12 with Eq.29, 
Eqs.35 - 37 yields an acceptable approximation. As is 
obvious from Fig. 4, than closer energy to the thresh- 
old of n particle production, the better approximation 
Eq.29. Therefore, if we choose the range of low energies, 
for example, up to 100 GeV, because in this range to- 
tal cross-section growth is observed, it is expected that 
the considered approximations will be acceptable for the 
large numbers of particles than those for which they were 
tested. In addition, as it follows from Fig.l0(b)-10(d), 
the accuracy of approximation Eq.37, as expected, in- 
creases with the growth of n. 
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FIG. 11. Comparison of the values of ^ G n \Pj ) obtained with consideration of all interference contributions (circles) and 

Z j =Z l 

the approximate values of (a n (zi)) p(zi-i) Az (crosses) for 11(a) - for n = 8 at y/s = 10 GeV, 11(b) - for n — 8 at y/s = 100 
GeV, 11(c) - for n = 9 at ^ = 10 GeV, 11(d) - for n = 9 at ^ = 100 GeV. 



Thus, within the framework of examined approxima- 
tions is possible to calculate the interference contribu- 
tions at sufficiently large n, and we can consider the de- 
pendence of total inelastic cross-section on energy y/s in 
the simplest case of multi-peripheral model taking into 
account all significant interference contributions. 

V. THE MODEL OF DEPENDENCE OF HADRON 
INELASTIC SCATTERING TOTAL CROSS-SECTION ON 
ENERGY 

Let us consider the magnitude 

ff ffi (^)="fLX(^ (39) 

n=l 

which within the framework of the discussed above model 
is an analogue of total inelastic scattering cross-section. 
Here n max is the maximum number of secondary parti- 
cles allowed by energy-momentum conservation law and 
L is the dimensionless coupling constant, which we con- 
sidered as a fitting parameter (see Eq.32 [2]). Since the 



calculation of a' n up to n = n max takes a long time, so 
in practice we restrict the upper bound of summation by 
those values of n, beyond which the neglected contribu- 
tions known to be smaller than the experimental error of 
cross-section measurements. 

The constant L can be fitted so that the dependence 
cr^ (y/s) looks like the behavior of total hadron-hadron 
scattering cross-section with a minimum about y/s = 10 
GeV. The result of such a fitting is shown in Fig. 13 (in 
that calculations we take proton mass as mass of primary 
particles and pion mass as mass of secondary particles). 

Quantitative comparison with experimental data re- 
quires the consideration of more realistic model than the 
self-interacting scalar (j) 3 field model. 



VI. CONCLUSIONS 

From obtained result, one might conclude that the con- 
sidered in [1] mechanism of virtuality reduction at the 
constrained maximum point of multi-peripheral scatter- 
ing amplitude may be responsible for proton-proton total 
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FIG. 12. The partial cross-section dependence on energy yfs calculated over all interference contributions (solid line) and by 
Eq.12 with the application of approximations Eqs.29, 35, 37 (dashed line): 12(a) - cr%(yfs)', 12(b) - a' 9 (^/s); 12(c) - 0io(y/s)] 
12(d) - 0n(i/s); 12(e) - (J\\(\fs). This approximation is acceptable at least in the range of parameters in which they are can 
be verified. 



cross-section growth when all the considerable interfer- 
ence contributions are taken into account. 

Just the revelation of mechanism of cross-section 
growth we consider as the main result of earlier papers [1 
and 2] and present work, since this mechanism is intrinsic 
not only to the diagrams of the "comb" type, but also to 
different modifications of considered model. 

Application the Laplace method allow to calculate an- 
other types of diagrams corresponding to various scenar- 
ios of hadron-hadron inelastic scattering and compare it 
with experimental data. 
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FIG. 13. Theoretical dependences of the a' 1 \yfs) 13(a) and cr /E (^/s) 13(c) obtained for the energy range y/s =1-7- 100 Gev 
at L = 5.51. First minimum for the total cross-section can be obtained only when we take into account contributions from 
the high multiplicities. Experimental data for the inelastic 13(b) and for the total 13(d) pp scattering cross-section [7 and 
8] presented for qualitative comparison with the prediction from our model. Note: data-points for the inelastic cross-section, 
obtained from the definition a ine i = (J total - & elastic- 



TABLE I. Results of numerical calculations of the eigenvalues of matrix D y . 
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